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^ . Abstract 

' In this talk we present some studies in the approach to equiHbrium of the classical Ai/)'' theory on 

■ the lattice, giving particular emphasis to its pedagogical usefulness in the context of classical statis- 
tical field theory (such as both the analytical and numerical evaluation of correlation functions) , as 
well as in the context of statistical mechanics (such as the equivalence of ensembles in the thermody- 
namic limit). Since the quartic term could be regarded as a smooth perturbation to the (integrable) 
gaussian theory, we also discuss the connection of our results with the KAM theorem, showing that 

■ the ergodic and non-ergodic regimes are not sharply separated. 

o 

^ ! 1 Introduction 



There has been some interest recently in classical fields as effective theories for finite temperature quan- 
tum field theory (see, for instance, |[l|, |2|). Furthermore, phase transitions, thermalization, ergodicity 
' and all that from statistical mechanics are still topics of much research in the context of field theories. 
This talk will be based mostly on Aarts' paper ||3|, which presents some results in the thermalization 
^ ' of the A(/)^ theory on the lattice for various ensembles. We will be interested, in particular, in deriving 
I analytical results related to the canonical ensemble in equilibrium and comparing them with the results 
of lattice simulations out-of-equilibrium, showing how the thermalization of the theory depends on A. 
This latter result is discussed in connection with the KAM theorem |^, |5|, which in the context of the 
Fermi-Pasta-Ulam problem has already been studied by Parisi elsewhere [||]. Although the simulations 
are done in a closed microcanonical system, in the thermodynamic limit we recover the equivalence 
of the ensembles @, therefore justifying the derivation of the correlation functions with the canonical 
probability measure. 
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2 The Classical A(/)^ Theory 

We will adopt the traditional A^^ action in 1+1 dimensions: 



/*/o 



dx 



(1) 



The equation of motion for this action could then be obtained directly using the field version of the 
Euler-Lagrange equation: 



dC d dC d dC 



d(p dx d{d^(j)) dt d{dt(p) 



= 0, 



where £ is the lagrangian density: 



\2 1.^,2/2 1\.4 



Then the equation of motion is: 



(2) 

(3) 
(4) 



The energy density is also easily obtained as the Legendre transform of C: 



whence we obtain the energy functional: 

£;[^,7r] = I dx 
Jo 



(5) 



with vr = dC/d(f) = dt4> the conjugate momentum to (p. 

3 Equilibrium Properties 

We start our statistical study by writing down the equilibrium partition function for the classical scalar 
field theory: 

Z = J dir J d(l)eicp{-PE[(l),Tr]), (6) 
where the functional measures dn and dcj) are the continuum limiting cases of, respectively, 

c/tt = JJ diTi, d(f) = Y\_ d(f)i 

i i 

Since the A<^^ term makes our <p integral non-gaussian, we shall analyze expectation values for the 
(gaussian) canonical momentum tt: 



(7r(x)), 



J dn J d<pex.p {—(3E[<p,7r]) 7r{x) 



2 



{7r(x)7r(y)). 



JdirJ dcpexp {-I3E[(I), n]) 7r(x)7r(y) 
Z 



which, after isolating the (p integral and canceling common terms with Z, reduce to the following 
expressions: 

_ JdiT exp(-/?£^o[7r])7r(x) / f ^ f -.s _ I direxp (-/?^o[7r]) 7r(x)7r(y) 

{7r[X))rp — - , [7T[X)TT[y))rp — - 



with: 



-^oM = / dx-n'^ and Zq = [ dvr exp (— /3i?oM) 
Jo 2 J 



In order to compare these expectation values with the numerical ones (i.e., in the lattice), we evaluate 
the above expressions using the discrete version with lattice spacing a: 



(vr(x)) 



1 

T - ^ 



d-Ki exp 



(7r(x)7r(y))^ = ^ 
Zq 



JJ^ diTi exp 



Let us begin with the expectation value {■k{x))j,, which after some simplifications with Zq reduces to: 



(vr(x)), 



S^!^ d-Kx exp i-^Trf) TTx 



j^^dl^xG^vi-^T^l 



0, 



(7) 



where we have used the fact that we are integrating an odd function in a symmetric interval. The 
two-point correlation function is obtained likewise: 



TTxTT. 



r r+oo 

J J—QO 


dTTxd-Ky exp 






TTxTTy 


II- 


diTxdTTy exp 









(8) 



where we took A;^ = 1 so that /3 = 1/T and have used the following general result for gaussian integrals: 



{VnVr. 



■y I- r+co 



2 X!^' 



\ / 1 



(9) 



with Aij a symmetric matrix and Z a normalization constant (see, for instance, Q). 
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As in 1^, we shall study observables constructed from the following two generating functions: 

G^'\t) = ^Y.^H^,t), G(^\t) = ^Y.^\x,t) (10) 

X X 

We now present a useful result for the correlation functions constructed from these generating func- 
tions in equilibrium: 

which could be easily derived using the following relation for gaussian integrals 0: 

{VnVmVpV,) = (a-') (A-') + (A-') (A-') + (A-') (A-') , 
\ / nm \ / pq \ / np \ / mq \ / nq \ / pm 

where A is the same matrix which appears in (|9|). 

We emphasize that these results were obtained in the canonical formalism, and are therefore restricted 
to either a system coupled to a heat bath, or an isolated system in the thermodynamical limit, where 
the canonical and the microcanonical ensembles are equivalent. Otherwise these results should be taken 
only as an approximation which becomes better for larger ensembles, i.e., close to the N ^ oo limit. 



4 The Lattice Formulation 

We have obtained the lattice version of the continuum Xcj)^ theory using second-order approximations 
for the derivatives in the equation of motion (^) (incidentally I also discuss these discretizations with 
some details elsewhere Q). We have then the leap-frog scheme: 



b^' = P (C-fi + C-i) + 2(1 - p)cP7 - - At^ 



(12) 



with p = (At/Ax)^. The above equation could then be solved iteratively for later times explicitly in 
terms of known values at times n. 

We have solved equation ([l2|) with periodic boundary conditions, together with the following (out- 
of-equilibrium) initial values for (p and vr (again from [^): 

(rie) 

(p{x, 0) = 0, 7r(x, 0) = ^Acos{27rkx/L - Vfc), (13) 

k 

where ipk is a random number in the range [0, 27r), and the sum over k has terms such that we only 
excite momenta of the order of the mass. The above initial conditions could be implemented using 
a first-order retarded difference scheme approximation for vr = dt(p, furnishing us with the following 
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discrete variables: 

(rie) 

</,o = 0, = -At^ Acos(27rA;x/L-'(/'yk)- (14) 

k 

We have also tried a number of different ratios At/ Ax in order to guarantee the stability of the 
solutions since the standard Courant condition ||8| cannot be naively applied here (due to the additional 
terms of first and third order). We have found stable solutions for p = 1/4 with m = 0.001 and A ranging 
from « to 1000. 



5 Observables 

We simulate our microcanonical ensemble by generating N^. "copies" of the system with a given fixed 
energy E, which depends only on the initial amplitude A, on the linear dimension L and on the number 
of excited modes rig, as can be seen by equations (||) and (13). 



In order to have a time-dependent expectation value for the observable O in this ensemble, we take 
the average over the copies at each instant t: 



n=l 



Since we are interested in the approach to equilibrium, we shall choose observables that could be 
easily compared with the equilibrium ones. In our simulations we have restricted ourselves to the function 
dev(7r) defined in |^, which is a measure of how much vr deviates from being gaussian. Using the already 



dev(7r) = ^ '-^ 1, (15) 



shown relation (11) and the following definition, we see that dev(7r) is a good measure of "how far" we 
are from equilibrium: 

with dev(7r) — > as the thermalization is attained. 

6 The KAM Theorem 

A useful theorem in the study of the ergodicity of finite systems is the KAM theorem (see, for instance, 
iS Hi)- We will follow these two references and stick to the qualitative view of the theorem, which 
roughly says that if the Hamiltonian of our finite system has an integrable term Hq (i.e., it has some 
integrals of motion) plus a smooth perturbation gH, 

H = no + gn, 

then for sufficiently small g the system will not be ergodic. This might suggest that there is a certain 
critical value of g which separates sharply the ergodic and non-ergodic regimes, and indeed it was found 
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Figure 1: Time-dependence of dev(7r) for different values of A. From top to bottom: A = 1000, 100, 10, 1. 
These results were run for an ensemble of Nc = 10, all of which with N = 200 lattice points. In order 
to make the graph clearer, we have reduced the fluctuations by taking local time averages of length 50 
time-steps. 

in some numerical simulations of the FPU problem Q. However, as we will show later, our system 
presents a smooth divergent behavior of the relaxation time as A — > 0, thus excluding the existence of 
the above critical value. 

The gaussian model, which is obtained by taking A = in our theory, has a closed-form solution with 
constants of motion and is therefore restricted to sweep a very limited area in the phase space. Thus, 
the ergodicity condition 

lim /* dt'g[C{t')] ^ 1 / dt,iC)5 {E - H{C))g{C) (16) 

*^00 JO J 

with dfi{C) the measure of the configuration space and O a normalization constant, is not satisfied, so 
this system is not a "good" statistical ensemble. 

However, in the non-linear A 7^ case (i.e., with a term of order (j)^ in addition to the linear term in 
the equation of motion), we might expect ergodicity for some values of E and A, and in fact this system 
has already been studied also by Parisi in the context of a chain of anharmonically coupled oscillators 
(the well-known Fermi-Pasta-Ulam problem ||6|). It is found in this latter paper that the system doesn't 
present the critical g said above; on the contrary, it always relax towards equilibrium, although with 
divergent time scales as y — > 0. Here we shall also analyze these issues, albeit with a different approach. 

7 Simulation Results and Discussions 

We have run the simulations on the lattice described previously for different values of A and (we have 
fixed both the lattice length L = 1 and the mass term m = 0.001, since these values together with A 



^,,.,^v^'#i(>JWA*^^^ ,^ 
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Figure 2: Behavior of the slope (3 (see Figure ^) for different values of A. We see here both the calculated 
values (circles) and the linear fitting. Since (3 measures the rapidity with which the thermalization is 
attained (at least initially), we see that the relaxation time diverges as A — > 0. 

ranging from 1 to 1000 allowed the system to behave in very different ways, thermalizing within t = 100 
for A f« 1000 and being almost locked in the initial phase for Awl). The results for different values of 
the coupling constant are shown in Figures || and || (the slope /3 used in Figure |2| is explained in Figure 
|3|). We notice that the relaxation times depend strongly on the coupling constant, being the faster times 
for greater A. 

These results are in good agreement with |^], which found no critical value for the coupling constant 
(albeit in the already mentioned context of the FPU problem) and also discuss the apparent existence of 
this critical value when one interprets naively the relaxation times in similar plots. We see from Figure |2| 
that the "rapidity" /3 vanishes in a power-law as A — > 0. Since in all of our numerical studies the function 
dev(7r) presented an almost monotonic increasing behavior, we found that /3, although not a precise 
measure of the relaxation time, has proven to be a very elucidative tool in estimating the time-scales in 
our restricted spectra of parameters. 

We have also studied our system with regard to its dependence on the number of sites N . We 
remember that dev(7r) was defined using canonical expectation values, whereas our system is a completely 
closed (microcanonical) ensemble. We justify this approach by noting that in the thermodynamic limit 
(i.e., N oo) the statistical ensembles should be equivalent (except, perhaps, for very rare systems which 
we won't touch here). We illustrate this behavior in Figure where we can see that our microcanonical 
dev(7r) tends to the canonical zero value as — cxd. The values used in this figure were obtained by 
doing an average in the "quasi-stationary" regime (following Aarts), where the system doesn't show any 
tendency of changing the expectation value of dev(7r). 
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Figure 3: How to obtain (5 from Figure ||. The symbols represent the numerical results (with local time 
averages), whereas the lines are linear fittings. The "rapidity" /3 is defined as the slope of the linear 
fitting from t = to 10. The upper curve is for A = 100 and the lower one is for A = 1 (the lower curve 
has been shifted down by 0.4 in order to make the graph clearer). 
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Figure 4: dev(7r) in the quasi-stationary regime (180 < t < 200, A = 1000) for different values of N . 
The triangles show the numerical results whereas the lines are straight-line interpolations. 
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8 Conclusions 



We have studied numerically the approach to equilibrium of the theory, showing the strong depen- 
dence of the relaxation times on the coupling constant A. In order to analyze how far the theory is from 
equilibrium, we have adopted an observable which could be defined locally in time during the simulation 
and also easily calculated in equilibrium using the standard canonical probability measure, depending 
essentially on the conjugate momentum vr. 

We have also defined a quantity (5 (which we called "rapidity") which measures how fast our theory 
is approaching to equilibrium. This quantity, as defined here, does not measure the overall behavior of 
the thermalization, but furnishes us with a good illustration of the divergent behavior of the relaxation 
time as the coupling constant vanishes. This result is not, at first sight, to be expected when one looks 
at the KAM theorem, since it might suggest a critical behavior of the relaxation times with A. 

We finish our study by showing how the system behaves for large values of the number of sites N , 
finding the vanishing of the quasi-stationary expectation value of dev(7r) as A?" — > oo. This latter result 
justifies the use of the canonical probability measure to calculate the equilibrium expectation values. 
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